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Abstract 

I discuss radiative corrections to the BFKL equation for high energy cross sec- 
tions in perturbative QCD. Due to the gluon Reggeization in the next-todeading 
Ins approximation, the form of the BFKL equation remains unchanged and the 
corrections to the BFKL kernel are expressed in terms of the twodoop contribution 
to the gluon Regge trajectory, the onedoop correction to the Reggeon-Reggeon- 
gluon vertex and the contributions from two-gluon and quark-antiquark produc- 
tion in Reggeon-Reggeon collisions. I present the results of the calculation of the 
BFKL kernel in the next-to leading logarithmic approximation, the estimate of the 
Pomeron shift and the next-to-leading contribution to the anomalous dimensions of 
twist-2 operators near j = 1. 
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1 Introduction 



The BFKL equation |l] is very popular now, mainly due to recent experimental results on 
deep inelastic scattering of electrons on protons obtained at HERA M , which show growth 
of the gluon density in the proton with decreasing of the fraction of the proton momentum 
carried by gluon. It can be used together with the DGLAP evolution equation || for the 
description of structure functions for the deep inelastic ep scattering at small values of the 
Bjorken variable x (see, for instance, 0] and references therein). The equation was derived 
for scattering amplitudes in QCD at high energies y/s and fixed momentum transfer \f—i 
in the leading logarithmic approximation (LLA) which means collection of all terms of 
the type [a s lns] n . This approximation leads to a sharp increase of cross sections with 
c.m.s. energy yfs. In fact, calculated in LLA, the total cross section a^ A grows at large 
c.m.s. energies as a power of s: 

-.LLA b /-p 
a tot T= . (1. 



/In s 

where Up is the LLA position of the most right singularity in the complex momentum 
plane of the t-channel partial wave with vacuum quantum numbers (Pomeron singularity), 
given by 

a 2 

= *-JVln2 (2) 

2 

for the gauge group SU(N) (N = 3 for QCD) with gauge coupling constant g (a s = |-). 
Therefore, the Froissart bound o to t < const(\as) 2 is violated in LLA. The reason of the 
violation is that the s-channel unitarity constraints for scattering amplitudes are not 
fulfilled in this approximation. The problem of unitarization of LLA results is extremely 
important from the theoretical point of view. It is concerned in a lot of papers (see, for 
example, || and references therein). 

The violation of the Froissart bound means that LLA can not be applied at asymp- 
totically large energies. But in the region of energies accessible for modern experiments 
it seems that the most important disadvantage of LLA is that neither the scale of s nor 
the argument of the running coupling constant a s are fixed in this approximation. These 
uncertainties diminish the predictive power of LLA, permitting to change strongly numer- 
ical results by changing the scales. From the practical point of view, since the results of 
LLA are applied to the small x phenomenology, it is extremely important to remove these 
uncertainties. Another important problem is the determination of the region of energies 
and momentum transfers where LLA could be applicable. To solve these problems we 
have to know radiative corrections to LLA. 

Therefore, the radiative corrections to the BFKL equation are very important, as they 
give the possibility to fix the argument of the running coupling, to define the scale of 
energy and to determine the region of applicability of the results obtained. My talk is 
devoted to these radiative corrections. 

The outline of the talk is the following. In Section 2, I remind the derivation of the 
BFKL equation in LLA and the solution of this equation. In Section 3, I discuss the 
general form of corrections in the next-to-leading logarithmic approximation (NLLA) and 
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present the two- loop correction to the gluon Regge trajectory and the contributions to the 
NLLA kernel from the one-loop correction to the one-gluon production and from the two- 
gluon and quark- ant iquark pair production in the Reggeon-Reggeon collisions. In Section 
4, all the contributions are collected together, the cancellation of infrared singularities 
is performed, the estimate of the shift of the Pomeron intercept and the next-to-leading 
contributions to anomalous dimensions of twist-2 operators near j = 1 are presented. 



2 The BFKL equation in LLA 

Despite the fact that the BFKL equation was obtained more than 20 years ago, till now 
a simple derivation of this equation does not exist, though attempts to do it continue 
(see, for examlpe, |§). In the original derivation Q the key role was played by the 
gluon Reggeization in QCD §]. In fact, the derivation can be performed without large 
difficulties if we adopt the Reggeization. It is worthwhile to stress here, that the gluon 
Reggeization in QCD means something more than usually assumed. Namely, it means 
not merely that there is the Reggeon with gluon quantum numbers, negative signature 
and trajectory 

j{t) = l + u(t) (3) 

passing through 1 at t — 0, but also that only this Reggeon gives the leading contribution 
in each order of the perturbation theory to amplitudes with the gluon quantum numbers 
in channels with fixed momentum transfers. Due to this property it is not difficult to 
calculate the leading contributions to the imaginary parts of elastic scattering amplitudes 
with arbitrary quantum numbers in t-channel at large s and fixed t using the unitarity 
condition. Full amplitudes are easily restored through their imaginary parts. The BFKL 
equation emerges from the representation of the amplitudes in the particular case of the 
forward scattering with vacuum quantum numbers in t - channel. 

For the elastic scattering process A + B — > A' + B' the gluon Reggeization means that 
at large s and fixed t, with 

s = (pa+Pb) 2 , t = q 2 , q = PA-PA>, (4) 
the amplitudes with gluon quantum numbers in t-channel have the factorized form 



i a \ A ' B ' -v i S 

K^SJAB - 1 A'AJ 



^B'B 5 (5) 



where T l A , A are the particle-particle-Reggeon (PPR) vertices. In LLA for the deviation of 
the gluon trajectory from 1 we have || 

where t = q 2 ~ q\ and D = 4 + 2e is the space-time dimension. A non-zero e is 
introduced to regularize Feynman integrals. The integration in Eq.(^) is performed over 
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(D — 2)-dimensional momenta orthogonal to the initial particle momentum plane. The 
PPR vertices can be presented as 



r AIA = g(A\T\A)T A , A , (7) 

where (t4'|T 1 |A) stands for a matrix element of the colour group generator in the corre- 
sponding representation (i.e. fundamental for quarks and adjoint for gluons). In LLA the 
helicities Xp of each of the scattered particles P are conserved, so that in the helicity basis 
we have 

^A'A = ^A'A = 5\ A ,\ A ■ (8) 



The s - channel unitarity relation for the imaginary part of the elastic scattering 
amplitude A\® can be presented as 

1 oo „ 

Im 8 Aii' = ^EE/ Kl +n d*AB + m (9) 

n=0 {/} 

where Im s stands for the s - channel imaginary part, A A u ^ +n is the amplitude of production 
of n particles with momenta ki, i = 1, ....n in the process A + B — > A + B + n, J2{f} means 
sum over the discrete quantum numbers of the final particles in this process, d^ A is 
the element of the final particle phase space. We admit all particles to have non zero 
masses (reserving the possibility to consider each of them as a group of particles) and use 
light-cone vectors p\ and p 2 such that momenta of the initial particles A and B are equal 
Pa = Pi + (^a/ s ) Pz and ps = P2 + ( m R,/ s )Pi respectively and s = 2(pip 2 ). Using the 
Sudakov decomposition 

h = Api + + k i± , saifii = k 2 - k 2 iL = k\ + k 2 , (10) 
where the vector sign denotes (here and below) transverse momenta, we obtain 

o 2 n+1 2 n+1 

S S i=0 S i=0 



n+1 , d/3 n+1 da A d fr tt d D ~ 2 k t± 

i=0 

with the denotations 



x^(gM=nin^, (id 



Pa = k 0i Pb = K+i- (12) 



In the unitarity condition (|9]) the contribution of order s, which we are interested 
in, is given by the region of limited (not growing with s) transverse momenta of pro- 
duced particles. Only this region is considered in the following. Large logarithms come 
from integration over longitudinal momenta of the produced particles. Therefore in LLA, 
where production of each additional particle must give the large logarithm (In s), they are 
produced in the multi-Regge kinematics (MRK). By definition, in this kinematics their 
transverse momenta are limited and their Sudakov variables CKj and fti , i = -7- n + 1, 
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are strongly ordered (in another words, the produced particles are strongly ordered in the 
rapidity space). Let us take, for definiteness, that 

a n+ i > a n > a n _i.... > a , A) > Ai > As— > An+i- ( 13 ) 
In this case the 5— functions in Eq. fllTD give us 

a n+1 ~ 1, fa ~ 1 (14) 

and therefore 

ao .fci (15) 
s s 

In MRK the squared invariant masses = (k{ + A;.,) 2 of any pair of produced particles 
i and j are large. In order to obtain the large logarithm from the integration over A 
for each produced particle in the phase space ([□]), the amplitudes in the r.h.s. of the 
unitarity relation (0) must not decrease with the growth of the invariant masses. It is 
possible only in the case where there are exchanges of vector particles (i.e. gluons) in all 
channels with momentum transfers % i — 1 -1- n + 1, whith 

i— 1 n+1 i— 1 

ft = pa - = -(pb - - Api - oa-m - J2 k i±> It - Qi± = -Qi 2 ■ ( 16 ) 

3=0 l=i 3=0 

Due to the gluon Reggeization the amplitudes of such processes in LLA have simple 
multi-Regge form: 



where is some energy scale, which is irrelevant in LLA, 

B Sl = {h-! + h) 2 ~ sA-iai = + A;'), U = g 2 ~ -g+ 2 , (18) 

Pi 

uj(t) and r^, p are the gluon Regge trajectory and the PPR vertices given by Eqs.(H) and 
©j © respectively; T^.^ftjft+i) are the effective vertices of production of particles Pi 
with momenta qi — qi + \ in collisions of Reggeons (i.e. Reggeized gluons) with momenta qi 
and — q i+ i and colour indices Cj and c i+1 correspondingly. In LLA all produced particles 
Pi must be gluons; therefore, the masses of produced particles are equal to zero. The 
Reggeon- Reggeon-gluon (RRG) vertex has the form [[I]] 



72m (ft, ft+i) = gT c % +i e;(k)C»(q l+u (19) 

where Tjh^ are the matrix elements of the SU(N) group generators in the adjoint rep- 
resentation, di is the colour index of the produced gluon, e(ki) its polarization vector and 
h — ft — ft+i its momentum; 

CXft +1 , ft) = -ft - ft+i + Pi(/- + 2^) - p 2 (f±i + 2^1). (20) 

K»Pl PlP2 ^iP2 PlP2 
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The amplitude A^?B> n entering the unitarity relation (^) can be obtained from Eq. ([H]) 
by the substitutions A — > A' , B — > B' , q^ — > q\ = q± — q, where q = Pa~ Pa 1 — <?±- 

Let us introduce the operators Vr for projection of two-gluon colour states in t— 
channel in the unitarity condition @ on the irreducible representations R of the colour 
group and use the decomposition 



T c% i+1 (T^ )* =^RMV R \c i+1 ^ +1 ). (21) 

R 



We'll be interested in the singlet (vacuum) and antisymmetrical octet (gluon) represen- 
tations. For the first of them 

<c lC ;iPoic J+1 c: +1 > = ^f;^ 1 (22) 



<c,c:ip 8 |c m c: +1 > = ^ c c +ic ' +i , (23) 



and for the second 

fac i c' i fac i+ ic' i+1 

~N~ 

so that one can easily find 

N 

c Q = N, c 8 = — . (24) 



Using the decomposition (|21|) we obtain from fllTf ) 



E7S I+1 (*,%+i)(7 c % +i (*,gm))* = E( c ^l^l c mc: +1 )2(27T) D - 1 /C^(g;,g m ;g), (25) 

Gi R 

where the sum is taken over colour and polarization states of the produced gluon and 
K,l R \q h q i+ i; q) = - 2 ^d-i ~Q,Qi~ 0) 



(2tt) d -i V (« - £+i) 2 

The decomposition ([H]) corresponds to the decomposition of the elastic scattering 
amplitude A^g' in (g): 

= E (^C , (2T) 

where Ar is the part of the scattering amplitude corresponding to the definite irreducible 
representation R of the colour group in t -channel. It is convenient to consider its partial 
wave f R (u, q)AB defined by 

f R (u>, q)iE' = r % (-) ~" Im. (A R )**' . (28) 

J so S \Sq/ 

The amplitude itself is expressed trough the partial wave as 

(A^J = ± t'~ -*L- ((=£)" - r (±)") /«(., ,")£?', (29) 
2-n Js-ioo sm(7ra;J \\ s J \s J ) 



where r is the signature and coincides with the symmetry of the representation R. For the 
gluon representation the Born contribution must be added into the r.h.s. of (29). The term 
with t takes into account the contribution to the amplitude from the u -channel imaginary 
part. Pay attention that the only antisymmetrical representation contributing to the 
decomposition fl2"T|) is the representation with the gluon quantum numbers. Therefore, 
only for this representation the contributions of s- and u- channel imaginary parts do not 
cancel each other. It means that in each order of perturbation theory the amplitudes with 
the gluon quantum numbers in -channels are dominant. 

Let us calculate the contribution f^ioOjq)^' into the partial wave coming from 
production of n gluons. Using Eqs.([TT|), (|I"3|)-(|T5l), ([18|), we have 



and from Eqs. (0), (§), and (p5|) we obtain 



Si 



k 2 k < 



x 



== S ^ II *^(* , 9) (31) 

QlQn+lJ » W / 



The index z/ here enumerates the states in the irreducible representation i?, so that 

(c i c' i \V R \c i+1 c' i+1 ) = J2(cic'i\V R \u)(u\V R \c i+1 c l+1 ) (32) 

V 

and 



4¥ = YXzb rfe 1 (H^|cn + i< +1 ). (33) 



BB \ x BB' 

B 

The sum here is taken over the discreet quantum numbers of the states A, B. 

For the singlet representation the index v takes only one value, so that we can omit it 
and have 

{cc ' im = vfe (34) 

whereas for the antisymmetrical octet (gluon) representation the index v coincide with 
gluon colour index and 

< CC '|P 8 |a> = (35) 
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where f a b c are the structure constants of the colour group. Therefore, we obtain 



iPp = ^0^=(P'\P)S, plXp , (36) 



where Cp is the value of the Casimir operator in corresponding representation, i.e. Cp 
x ~~ 1 for quarks and Cp = N for gluons; 



2N 



IP? = -ig 2 ^-{P'\T a \P)5, lXp = -ig^-Y a pip . (37) 



N „. „ \/N 



The integration over Sj in Eq. fl3T|) is performed from some fixed (independent from s) value 
to infinity. Note, that the essential integration region in Eq. ( p9|) is u ~ (Ins) -1 , so that 
in LLA we can omit terms of the type uohisp, ulnk^, as well as Lo{tj) In sp, c<j(t-) lns#. 
It corresponds to the statement that in LLA the scale of energy is not fixed. Therefore, 
independently from the lower limit of the integrations over Sj we obtain 



n+l AD-2 



f(»V, -*A'B' _ 1 / TT » g'-L 



(2vr)^ J 11 (w _ „ (f .) _ W (tj)) 5 a(gj - g)2 



v \i=l 

Let us present the partial wave in the form 



M».1>& = qJf* I tf ( >& V &G ^' fa aigg> • (39) 

Then for the function G^' (which can be called Green function for scattering of two 
Reggeized gluons) we obtain 



Gl R \q A ,q B ;q) = j: f II 

n=0 J \i=l 



n+l ^7D-2 



(f[ q i+1] q^j q^(q A ~ q) 2 q 2 (q B ~ q) 2 5 D ~ 2 (q 1L - q A± )5 D - 2 (q n+1± - q B± ). (40) 

It is easy to see, that the sum in the r.h.s. is the perturbative solution of the equation 

wGg^ft, ft; q) = q x 2 {qx ~ &SP- 2 \& - S) 

+ / ^^- ^(a, g/; flG^tfi', g 2 ; g), (41) 

J q\ (?i -0 

where the kernel 

/C (R) (gi, g 2 ; g) = + w(( gi - g)i))gl 2 (gl - flW^gi - g 2 ) + /C^(gl, g 2 ; g) (42) 



consists of two parts; the first of them, so called virtual part, is expressed in terms of the 
gluon Regge trajectory and the second, related with the real particle production, is given 
byEq.(B§). 



Let us consider the partial wave f B \uj; q\, q) B for the Reggeized gluon scattering off 
the particle B defined as 

ff\u; £, q)% = I d D -*q BA .GWfr, q B ; $ *?*\ 2 , (43) 

J qB \qB - q) 

where the impact factor Igrg is given by Eq. (|33|) . This partial wave satisfies the equation 
uft\u; q u q) B B = igg + J ^^y ^^, & q)f { »\u; & q) B B . (44) 

It is easy to see that for the case of the gluon representation the solution of this equation 
is 

s ^^) b b =zMky (45) 

Therefore Eqs. ©, ©-(H) give us 

f (. , rf\A'B' _ r(8,a) 1 f d°~ 2 qi± r (8,a) 1 



(tor)"-* J tfiqi-q) 2 B ' B uj-uj(ql) 

- iA ' A t{u-u{t)) iB ' B - m 

Exactly the same result one gets substituting Eq.fl5|) in Eq. (|28|) . So, we have a kind of 
"bootstrap": we input the existence of the Reggeized gluon and obtained it as the solution 
of the equation derived from the unitarity. 

Strictly speaking, this "bootstrap", although being very impressive and meaningful, 
can not be considered as a rigorous proof of the gluon Reggeization. For such proof we 
have to reproduce not only the form (Eft of the elastic amplitudes, but also the Reggeized 
form ([17]) for all inelastic amplitudes. It was actually done in LLA in Ref. ||, so that in 
this approximation the proof of the gluon Reggeization really exists. 

Let us turn now to the most interesting case - vacuum quantum numbers in t -channel. 
From now, only this case will be considered, with additional simplifying restriction of the 
forward scattering, i.e. A' = A, B' = B, q = 0. Instead of the imaginary part of the 
forward scattering amplitude, we'll consider the total cross section ctab(s), 

(?ab{s) = • (47) 

s 

From Eqs.(|29|), (j39l) we obtain 

rS+ioo du 1 [ , D _ 2 jD-2 ( s 

aMs) = Leo 2ri (2^F=J J d ^ ^ U 
q A q B 
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where 

$ a(u) = if A = ^ N l_ x E |ry 2 ; u = -pa, 

®B(q B ) = I { BB= J-_ El r ^| 2 ; QB = -P B (49) 

^ N ~ 1 B,c 

and 

G u {q 1} ga) = ^2 \B (50) 

The PPR vertices Vp tp are defined in Eqs.(0), (§) and the Green functions Gffl{q\, q 2 \ q) 
in Eqs.(PD, (pE2|). Let us note that for quarks and gluons the impact factors $p(qp) 
don't depend really on qp in LLA. We indicate this dependence having in mind possible 
generalizations to the cases of scattering of other objects and corrections to LLA. The 
Green function Gw(q\, q 2 ) satisfies the equation 

00 GM, Q2) = S D ~ 2 (qi - &) + / d D ' 2 q K{q u q) G u {g q 2 ) (51) 

where 

&) = 0) = 2 Loi-q*) ^{q, - q 2 ) + K r {q %1 q 2 ) . (52) 

Here u(— q 2 ) is the gluon Regge trajectory given by Eq.flf^) and the integral kernel 
)C r (qi, q 2 ), related with the real particle production, is given by Eq. (|26|) at q = 0: 

r _ ££ 0) (<fi>92;0) 1 1 Gw m2 

Mft.fc) " ee ~ 2(27r) n- 1{N2 _ ci E Gi l7 clC2 fe, ?2 )l 



(2tt)^ (gW 2 ) 2 ' 



(53) 



Taken separately, the virtual (|6|) and real ([53]) contributions to the kernel (|52|) lead to 
infrared singularities. Indeed, 

whereas the real contribution fl53|) gives the term of order 1/e after integration around the 
point qi — q 2 = 0. These singularities cancel each other in Eq.(^TJ). Let us demonstrate 
it for the kernel averaged over the azimuthal angle between qi and q 2 . In fact, only the 
averaged kernel is relevant until we don't consider spin correlations because in this case the 
impact factors entering Eq.(|48|) depend only on squared transverse momenta. It's clear 
therefore that the high energy behaviour of cross sections is determined by the averaged 
kernel. 



Performing expansion in e and keeping only the terms giving non- vanishing contribu- 
tions at e — > in (pi), we obtain for the averaged kernel 



Instead of the dimensional regularization we can introduce the cut off \q 2 — q 2 2 \ > A 2 in 
the kernel /C r (<fi, $2) changing correspondingly the virtual part. We obtain 



g^N 

471-3 







-A 2 ) 




-> 2 -* 2 

Qi ~ Q2 





(56) 



The representation fl56|) permits to find easily such form of the kernel for which the 
cancellation of the singularities is evident. For this purpose it is enough to present 



2 In 



A 2 



dq 2 







■A 2 ) 




Qi 2 ~ Q 2 2 





(57) 



with 0(1) = 1. Evidently, the mentioned form is not unique. The one adopted in |l] and 
used in literature is 



Jd D - 2 q 2 JC(q u q 2 )f(q 2 2 



71 



-> 2 

2^1 , 



<? 2 2 l<? 2 2 



(£ 2 ) 2 + 4(g 2 2 ) 2 



<l2 ~ ?1 

) 



In the following we'll use another choice: 



J d D 2 q 2 JC(q u q 2 )f(q 2 



7T 



dq 2 



W2 2 - Qi 2 \ 



(5f 



(59) 



Of course, the representations (|58[) and (p9| ) are equivalent. Both of them make explicit 
the scale invariance of the kernel, due to which its eigenfunctions are powers of q 2 . We'll 
take them as (q^ 2 ) 7 " 1 an d denote the corresponding eigenvalues as ^f JL X B (l)'- 

f d D ~ 2 q 2 JC(q u q 2 ){q 2 V- 1 = ^*( 7 )(£»)?-\ (60) 

J 7T 



so that 

x B {i) = 2^(1) - ^(7) - - 7) , ^(7) = r'( 7 )/r( 7 ) . (6i) 

The set of functions (^* 2 2 ) 7 1 with 7 = 1/2 —00 < f < 00 is complete, so that we 

have 



s + io ° did 

8—ioo 2711 



+00 



dv 



00 27P 



X 



(62) 



The cross section exists only if the impact factors possess a good infrared behaviour; 
otherwise it turns into infinity. In fact, it is infinite for scattering of colour particles, as 
it should be, and finite for colourless ones, because for them 



<S>p{qi 



1p 



(63) 



for small qp. The maximal value of x B {l) 011 t ne integration contour in Eq . (|62"D is x(l/2) = 
4 In 2, that corresponds to the maximal eigenvalue of the kernel ujp = 4N(a s /ir) In 2. At 
uj = Up the partial wave has the branch point that leads to the growth ([!]) of the cross 
section. 
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3 The Next-to-Leading Approximation 



The next-to-leading logarithmic approximation (NLLA) means that all terms of the type 
a[a s ln(l/x)] n have to be collected. It was argued in Ref. WU that in this approximation 



we can use the approach which coincides in the main features with that used in LLA. 
In general, the programme of the calculations is analogous to that in LLA. The final 
goal is the elastic scattering amplitude. It has to be restored from its s- and u-channel 
imaginary parts. The s- channel imaginary part is given by the unitarity relation (|9|). 
Evidently, Eqs.(p^)-p9|) expressing the elastic scattering amplitudes in terms of their s- 
channel imaginary parts remain unchanged. In the multi-Regge kinematics (MRK) in 
NLLA, as well as in LLA, only the amplitudes with the gluon quantum numbers in the 
channels with momentum transfers do contribute. It was mentioned after Eq. (|29| ) that 
in each order of perturbation theory these amplitudes are dominant, because only for 
them there is no cancellation between s- and u - channel contributions. Moreover, for 
the same reason only in these amplitudes the leading terms are real, whereas in other 
amplitudes they are imaginary. Therefore, the appearance in the r.h.s. of the unitarity 
relation (H) of amplitudes with quantum numbers in tj-channels different from the gluon 
ones leads to loss of at least two large logarithms and therefore can be ignored in NLLA. 
Evidently, it is a peculiar property of NLLA. In the approximations next to NLLA such 
amplitudes do contribute. 

As before, the key point in the calculation of the amplitudes contributing in the uni- 
tarity relation is the gluon Reggeization. In MRK the real parts of the contributing 
amplitudes (only these parts are relevant in NLLA because the LLA amplitudes are real) 
are presented in the same form ( p9|) as in LLA. So, in this kinematics the problem is 
reduced to the calculation of the two-loop contribution u^ 2 \t) to the gluon Regge trajec- 
tory uj(t) and the corrections to the real parts of the PPR- and RRG- vertices. Let us 
note here, that the PPR-vertices, as well as in LLA, enter the expressions for the impact 
factors (see Eq. (P3|)), but not the expression for the kernel (see Eq.p^),([25|)), so that the 
corrections to these vertices appear at the intermediate steps of the calculations only. The 



one-loop corrections to the PPR vertices were calculated in Refs. [[11], 0. Though they 
are necessary for the calculation of the corrections to trajectory and RRG vertex, they 
don't enter the kernel explicitly, therefore I don't present them here. 

But in NLLA, contrary to LLA, MRK is not a single kinematics that contributes 
in the unitarity relation @. Since we have a possibility to loose one large logarithm 
(in comparison with LLA), the limitation of the strong ordering fllTf ) in the rapidity 
space can not be implied more. Any (but only one) pair of the produced particles can 
have a fixed (not increasing with s) invariant mass, i.e. components of this pair can have 
rapidities of the same order. This kinematics was called [|IIJ] quasi-multi- Regge kinematics 



(QMRK). We can treat this kinematics including, together with the one-gluon production, 
production of more complicated states in the Reggeon-Reggeon (RR) collisions, namely 
gluon-gluon (GG) and quark- ant iquark QQ sates, as well as production of "excited" states, 
i.e. sates with larger number of particles, in the Reggeon-particle (RP) collisions in the 
fragmentation region of one of the initial particles. Therefore, the partial wave (p8|) can be 
presented in the same form fl39|), but with modified impact factors and RR Green function. 
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In the definition of the impact factors ([J3|) we have to include radiative corrections to the 
PPR vertices and the contribution of the "excited" states in the fragmentation region. The 
equation (f4T|) for the RR Green function remains unchanged, as well as the representation 
( f42"D of the kernel, but the gluon trajectory has to be taken in the two-loop approximation: 

u(t) =io^(t)+uj^(t), (64) 

and the part, related with the real particle production, must contain, together with the 
contribution from the one-gluon production in the RR collisions, contributions from the 
two-gluon and quark- ant iquark productions. The one-gluon contribution must be calcu- 
lated with the one-loop accuracy, whereas the two-gluon and quark- ant iquark contribu- 
tions have to be taken in the Born approximation. In the following we consider only the 
case of the forward scattering and present the part of the kernel related with the real 
production in the form 

KM, = ICZg° OP (Qu <f 2 ) + K%£§ofa, <f 2 ) + Cqq(?i. (65) 

Here there is a subtle point. Calculating the two-gluon production contribution to 
the kernel (as well as the contribution to the impact factor from the gluon production 
in the fragmentation region) we meet divergencies of integrals over invariant masses of 
the produced particles at upper limits (let us call such divergencies ultraviolet). These 
divergencies correspond to the uncertainties of the lower limits of integrations over Sj (see 
Eq. (|31])) in MRK, which were not important in LLA, but can not be ignored in NLLA. 
Of course, the reason for the divergencies is the absence of a natural bound between MRK 
and QMRK. In order to give a precise meaning to the corresponding contributions and 
to treat them carefully we introduce an artificial bound (which, of course, disappears in 
final results). The discussion of the separation of the MRK and QMRK contributions 



presented in subsection 3.3 is based on the paper [[13 



In the next subsections we'll discuss various contributions to the kernel. 



3.1 The Reggeized Gluon Trajectory 

The two-loop corrections to the gluon trajectory were obtained in Refs. [O, [L5I]. They 



were expressed in terms of discontinuities of QCD scattering amplitudes with gluon quan- 
tum numbers in t-channel calculated at large energy y/s and fixed momentum transfer 
\/—t in the two-loop approximation. The processes of the quark-quark, gluon-gluon and 
quark-gluon elastic scattering were considered. Independently from the process, it was 
obtained for the case of nf massless quark flavours: 



where t = q 2 = —q 2 and 

m , = / 4^1* [.„ ( *LJ\ - 2^D - 3) 

A(2ti) u 1 J q 2 2 {q2-q) 2 [ \{qi - ? 2 ) 2 J 
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3-£ 
2 



+ 2 "<2 



+ 



2g 2 Nn f T (2 



(D-3)(D 
r2 (f) 



+ 



9 



4) 4(D-l)(D-3) 



(f-2) 



• . , (67) 
(4vr)fr(D) ' 

Eqs.(^6j) and (|6T| ) give us a closed expression for the two- loop correction to the gluon 
trajectory. The independence from the properties of the scattered particles, which appears 
as the result of remarkable cancellations among various terms, sets up a stringent test of 
the gluon Reggeization beyond the leading logarithmic approximation. 

The two-loop correction to the trajectory contains both ultraviolet and infrared di- 
vergencies. The former ones can be easily removed by the charge renormalization in the 
total expression for the trajectory. Since the trajectory itself must not be renormalized, 
we have only to use the renormalized coupling constant instead of the bare one g . In 
the MS scheme 



9 = 9^ 



1 + 



11 

y 



2_nf_ 
3 N 



2e 



where 



9lNT(l-e) 



( 47r )2+e ' 

and calculation of the integrals in (|66|) , (|67|) gives |H| : 



(68) 



(69) 



"(t) = -9l 




7T 



C(3) 6 2 - 



(70) 



The remarkable fact which occurred is the cancellation of the third order poles in e existing 
in separate contributions to the gluon part of (pEp. This cancellation is very important 
for the absence of infrared divergences in the corrections to the BFKL equation. As result 
of this cancellation, the gluon and quark contributions to u^it) have similar infrared 
behaviour. Moreover, the coefficient of the leading singularity in e is proportional to the 
coefficient of the one-loop (3 function. This means that infrared divergences are strongly 
correlated with ultraviolet ones. The correlation is unique in the sense that it provides 
the independence of singular contributions to u(t) from q 2 . Indeed, expanding Eq. ([70|) 
we have 




9u 



3 3iVVU 2 





404 _ , 56 n f 



(71) 
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Eq. ( |TTD exhibits explicitly all singularities of the trajectory in the two-loop approximation 
and gives its finite part in the limit e — > 0. Let us stress that the independence of the 
singular contributions to the trajectory from q 2 is necessary for the cancellation of the 
infrared divergencies in the BFKL equation. 



3.2 One-Gluon Production Contribution 



The simplest process for extracting the RRG vertex is production of one gluon with 
momentum k\ in MRK: 

s > Si,2, |*i,2| , s i = (h+PA) 2 i s 2 = (h+pz) 2 , sis 2 ~ sk 2 , t i = q 2 ~-q i 2 . (72) 



It is worthwhile to remind, that the MRK amplitudes beyond the LLA have a com- 
plicated analytical structure. For the one-particle production amplitude, assuming the 
Regge behaviour in the sub-channels s\ and s 2 , from general requirements of analiticity, 
unitarity and crossing symmetry one has (see Refs. 11, 16]) 



aAB+1 
AB 



-r 



Cl 



-Y c ? 



4~AA tl c 1C2h BB 



X 



Q \ U)\—U>2 

I — 

sr 



a \ &2 / q \ ^2 

SrJ \Sr 



R 



+ 



-s 2 

sr 



— S 



Sr 



Sr 



(73) 



where Ui = u(ti) and the right and left RRG vertices R^ SR ^ and L^ SrS) are real in all physical 
channels. Fortunately, as it was explained at the begining of this section, in NLLA only 
real parts of the production amplitudes do contribute , because only these parts interfere 
with the LLA amplitudes, which are real. Therefore, for our purposes we can neglect the 
imaginary parts and present the amplitudes in the form ( |T7D with 



L) 1- 




2 Sr 

where u^(t) is the one-loop contribution to the gluon Regge trajectory. 



(74) 



The right and left RRG vertices R and L were calculated in Refs. ||11|| . The calculations 
were performed in the space-time dimension D 7^ 4, but terms vanishing at D — > 4 
were omitted in the final expressions. Unfortunately, such terms should be kept, because 
integration over transverse momenta of the produced gluon leads to divergency at k\± = 
for the case D = 4. Therefore, in the region ki± — > we need to know the production 
amplitude for arbitrary e. The corresponding calculations were performed in |T7[ . 
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Let us note here, that whereas the dependence of the Regge factors from the energy 
scale sp was beyond the accuracy of LLA, in NLLA it has to be taken into account, 
though we can use an arbitrary scale. It means, that the RRP-vertices, as well as the 
PPR-vertices become dependent from the energy scale sp. In the following we'll show 
explicitly this dependence, denoting the PPG-vertices by and the RRG-vertices 

by 7ci?2 Gl (<?i> ft)- From the physical requirement of the independence of the production 
amplitudes from the energy scale we have: 

r<->" = r<'«»" (|) K r«> (1 + l r (t pP ) h (g) ) , (75) 

where tpp> is the squared momentum transfer from P to P', and 

* 7g? 2 )Gl (ft,ft) (l + ^("1 + ^)ln (^Yj (76) 



In Refs. JTT|,|jl3 the PPR vertices were extracted from the elastic scattering ampli- 
tudes assuming the representation for them; therefore, the scales —t P p were used for 
these vertices. In turn, the production vertices R and L were extracted from the one- 
gluon production amplitude using the representation (|73|), with the scale in the Regge 
factors equal the renormalization scale /1 2 , but with the PPR vertices taken for the scales 
—t P p. So, several scales were mixed in the one-gluon production amplitude. Though such 
mixing of scales is not forbidden, it is quite inconvenient when we consider production 
of arbitrary number of gluons. Therefore we'll use for definition of all vertices the rep- 
resentation ( [T7| ) with the single scale sp. Then, according to Eqs.(|75|), ([H]), the vertices 
7cic 2 ^ Gl (9i; ft) differ from given by Eq.flTiD with sp = /j 2 and vertices R and L taken from 



Refs. |TT|,[|l7j by their LLA values (|T9| ) multiplied by the factor 



~ U (1 >(*i)ln 



(77) 



Therefore we have 
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As for the vertices R( SR ^ and LS 8R ' separately, they can be found taking into account that 
the combination R — L does not contribute in the LLA (see ([T4|)) and therefore can be 
taken in the first nonvanishing approximation. We obtain 



R (sr) _ L (sr) =R_L, 



l + ^Mln 



In 



-t 2 



where R and L are taken from Refs. [pi , ||17|| . So, we have 
R (s R ) + L (s R ) =2g ^- %l{kl) (C(q 2 , qi ) ll + g 



11 Uf s 

~Q ~ 3N' 



(79) 
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One can check that substitution of (|80|), (|81| ) into (|74|) gives (|78|). 



*0) 
31) 



It will be shown in the next subsection, that for the appropriate separation between 
MRK and QMRK the contribution K-°^Q loop (qi, q 2 ) to the kernel is equal to 



one— loop 
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From Eqs. 



we obtain [B 
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3.3 Separation of the MRK and QMRK contributions 

Let us rewrite Eq. fl28|) for the case of the forward scattering: 

r°° ds f s\~" , . , iS 

f(u,s ) AB = — — cr AB (s), (84) 



using relation fliTD and introducing so as the argument of the partial wave because in 
NLLA the partial wave becomes dependent on the energy scale. For the contribution 
fMRK^t s o)ab into the partial wave coming from production of n gluons in MRK, re- 
peating the steps leading to Eq. ([H]) we obtain 

*o)ab = / ill <l l '-\±— (-) \ ~== I - " s a ) 




x 



/^1= ^ 11^0(^,^+1;^) — 2 . (85) 

V?n+l/ ?1 Vi=l / ?n+l 

In order to give the precise meaning to MRK we introduced here the parameter s\. 
By definition, MRK means that all squared invariant masses are larger than this fixed 
(not growing with s), but as large as it would be desired parameter. The contribution 
fMRK( u i s o)ab depends on this parameter although we don't indicate this dependence 
explicitly. The dependence from this parameter disappears in the final answer for the 



partial wave. The impact factors (c.f. (49)) 



$^(g P ; s R ) = -jJ== £ |r<£ )c | 2 ; QP = ~VP (86) 

acquired the upper index B + v for the denotation that they represent the sum of the 
Born (LLA) contribution and the virtual correction, as well as the additional argument to 
show their dependence from the scale s R used in the Regge factors in Eq . (|17D . Remind, 
that besides the presented contributions, the total NLA impact factors contain the con- 
tributions from the gluon production in the corresponding fragmentation regions as well. 
The part of the kernel 

ICrrgH, q i+ i] s R ) = 2 (2n) D - 1 (N*-l)q l 2 q l l 1 ^ G l7 Kr^' *+0| 2 - ( 8T ) 
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represents the contribution of the gluon production with the one-loop accuracy and also 
depends on sr. 

It is convenient to change the argument Sr for k 2 in /Crrg(<&> qi+i] sr) and for So in 
$p +v (qp; sr) using Eqs. (|75|) , (|76|) . Then after integration over s, we obtain: 
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In 



In the last equality we performed the expansion in (u — 2uj) and have taken into account 
the first terms of the expansion only, as it should be done in NLLA. Let us note for 
defineteness that for n — 1 in the last factor in Eq.(|55|) only the two last terms remain, 
whereas for n = the whole factor is equal unity, because in this case in Eq. (^) instead 
of integrations over partial Si we have an integration over the total s which has to be 
performed from sq to oo. 

Now let us turn to QMRK. Pay attention, that since the contribution of this kinematics 
is subleading, we can ignore the dependence from the energy scale. In the exact analogy 
with the PPR and RRG vertices let us introduce the effective vertices 7^+i( ( 7«> an d 
7^2+i(<?«5 Qi+i) f° r two-gluon and quark-antiquark productions in the Reggeon-Reggeon 
collisions as well as the vertices Tp, P for production of the "excited" state P* in the 
fragmentation region of the particle P in the process of scattering of this particle off the 
Reggeon, so that the amplitudes of production of n + 1 particles in the QMRK are given 
by Eq. ( |TTD with one of the "old" vertices changed for corresponding "new" vertex. With 
this definition the contribution fj^MRK (^0 ab of the n + 1 particle production in QMRK 
is: 

1 r f n +} d D - 2 < 
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where $^(<fp) and K^\q^ qi+i) denote the Born (LLA) values for the impact factors 
( f49"D and the part of the LLA kernel related with the real gluon production (|53"D, & P \qp) 
means the corrections to the impact factors due to the production of the "excited" states 
in the fragmentation regions and 1CQMRK(<lj,<lj+i) appears as the contribution to the 
kernel from QMRK. The last two values depend on the boundary s\ between MRK and 
QMRK, though we don't indicate this dependence. Taking into account that the phase 
space flu]) includes d D - 1 fc i /((27r) D - 1 2e, ; ) and that d D h = dk 2 d D - 1 k i /(2e i ), we have 

£>QMRK(qj,qj+i) = 2 (jV2 _ \)q 2 q 2 +1 / d P p ^ e ( Sk ~ 

X E 17^(^+1) I 2 - (90) 

Cj,Cj + l,Pl,P 2 

Here the sum is taken over colours of Reggeons and all quantum numbers of produced 
particles. The produced particles P\P 2 can be quark-antiquark as well as two gluons; 
denoting their momenta by l\ and l 2 we have 

2 d £) - 1 l i 

dp PlP2 = dk)b{kj l 2 ) \{ ( 27r ) D -i2 e . - ( 91 ) 

Note that in the case of the two-gluon production the gluon identity must be taken into 
account by the limitation on the integration region. We can represent Eq. ([X]) as the 
following: 

£QMRK{qj,q j+ i) = {N - 1) / ds RR (0 , Do -? 2 - r2 > ( 92 ) 

\ Z7T ) z qj qj+i 

where (TR R _ 2 {k 2 ) is the total cross section of the two-particle production in collision of 
Reggeons with momenta qj and —qj+i, (qj — <?j+i) 2 = srr, averaged over colours of the 
Reggeons, and Irr is the invariant flux: 



irr = V (srr - q] - q] +1 ) 2 - ^qjqj+i = 1\j (Wj+i) 2 ~ QjQj+i; q\ = ~& '■ (93) 

Analogously, the correction $ P \qR.) to the impact factor due to the production of the 
"excited" states in the fragmentation region of the particle P at collision of this particle 
with momentum p P and Reggeon with momentum — q R , {p P — q R ) 2 = s PP , is given by 



4 r) (fe) = E / dp G p9(sA - s PR )\T 



GAA\ 



= jw^i I ^ 2Wp ^ p ; (g 7 )g(gA ~ gpfl) . (94) 

J f 2tt)2s 

Here we have taken into account that, with the normalization used the matrix element 
for the RR — > P1P2 transition coincides with the corresponding effective vertex, whereas 
for PR — > P* differs from such vertex by factor \/2s. 
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Note, that with the NLLA accuracy the corrections obtained by the expansion in 
{u — 2u>i) in Eq. d88|) for the case of the (n + l)-particle production can be presented in 
the same form as the r.h.s. of Eq. (|89|) with the substitutions: 

£^ _ _ / ^5^®,, |)i in ( , 
<1r J q 1 \{qR-q) 2 s j 



Kqmrk%, Qj+i) ~^ ~ J d ° 2 v±fci B Xqjiv)fci B) (tqj+i)2 ln 



A 



(95) 

Therefore from Eqs. flS8|)- (|90| ) we obtain that with the NLLA accuracy the total partial 
wave can be presented as 



where K, r (qi,q 2 ) has the form (|65| ) with JC° R n ^ G loop (qi, q 2 ) given by Eq. (|82|) , 

^tt)^ 2 ^ 

-Born ^ ^ \ _ / AT 2 -, \ f j„ %1 ' RRO RR^Gg(s Rr)9(s \ — S RR 

5 RR 



(96) 



- / dP-^JC^i&tiK? \lq j+ i)\\n ( g-) , (9f 

J 2 \{Qj - q) 2 {qj+i- q) J 

and the impact factors have the representation 

(2ir) D -" P " f 2 (,-H-!) 2 V(&-«%, 
Here we used the form (|53| ) of /C r B )(gi, ^). 



(99) 



In the last two equations the MRK boundary sa have to be taken much larger than 
typical squared transverse momenta, so that the dependence from s\ disappears due to 
the factorization properties of the Reggeon vertices in the regions of strongly ordered 
rapidities of produced particles. Namely, the two-gluon production vertex is given by the 
product of the RRG vertices: 

Tg&fo, Qi+i) - 7§(fc, Qi ~ h) ^jjT§ +1 (?< - h, q j+ i) (100) 

at (pBh) <C (pbIi), (paIi) (pa^), where Zi and Z2 are momenta of the produced gluons, 
and the vertices for the gluon production in the fragmentation regions are given by the 
product of the corresponding vertices without gluon and the RRG vertices: 

t gaa * T AA-^jy^c(qR + 1, q R ) (101) 
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at (p B l) < (pbPa), (PaPa) < (paOj ^ is tne gluon momentum, q R = p A - p A - I ; 

1 



(102) 



at (p A l) < (paPb), (PbPb) < (PbO , <1r = Pb ~ Pb + I ■ 



3.4 Two-Particle Production Contributions 



Investigation of the two-gluon production contribution was started in Ref.|T0|; the next 
step was done in Ref.|18|. The final result was obtained in Ref.[jl9| by calculation of 
differential cross section of the two-gluon production in the RR collisions and integration 
of this cross section over relative transverse and longitudinal momenta of the produced 
gluons. In a suitable for integration form we obtain: 
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(103) 



where l x and Z2 are momenta of gluons produced in collision of Reggeons with momenta 
q x and -q 2 , l x + l 2 = q x -q 2 = k x , 

qi = P1P1 + qi±, <?2 = -a 2 P2 + q 2 ±, U = q 2 ± = -q 2 , 



I 2 

h = x(3 x p x + — 5— p 2 + h±, 
xp x s 

tx = (qi - h) 2 = — ((h - xq x ) 2 + x(l - x)q{ 
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(h + k 



{(l-x)h-xl 2 ) 



-(1 - x)l* 



xl 2 ■ 



(104) 



x(l — x) 

Here we omitted terms which vanish after integtation. The intermediate parameter 5r is 
introduced in ( |103| ) to exclude divergencies of separate terms. This parameter should be 
considered as infinitely small. In the total expression for the integral the dependence on 
5r vanishes. 

Note, that taken separately some terms in the r.h.s. of Eq. (|103|) contain "ultraviolet" 
divergencies at large but these divergencies are artificial in the sense that they cancel 
each other, so that the total integral is convergent at large|Zi|. In the infrared region for 
e — > all terms have no more than logarithmic divergency. All divergencies, "ultraviolet" 
and infrared, can be regularized by the same e. Since in the BFKL equation the kernel 
has to be integrated over q 2 , or equivalently, over ki, and the kernel is singular at ki = 0, 
we should be careful in performing the expansion in e , in order to keep all terms which 
give nonvanishing contribution in the physical limit e — > after integration over k\. The 
result of the integration in ( 103 ) is: 
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The r.h.s. of Eq.( [105| ) exhibits several terms with unphysical singularities, such as 
at q 2 = q 2 and (q\ + <f 2 ) 2 = 0. It is not difficult to see, that the only real singularity 
is at A; 1 2 = 0. All others are spurious and singular terms cancel each other. As for the 
singularity at k^ = 0, Eq. ( |105| ) contains all terms which give in the physical case e — ► 
non-zero contributions after the subsequent integration over k\. For the terms which are 
infrared divergent at D = 4, the region of such small ki, that is ln(l/fc 1 2 ) ~ 1/e , gives 
essential contribution in the integral over k\. Therefore, we can not expand (fc 2 //i 2 ) e in 
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powers of e in such terms. Note, that the calculated contribution is explicitly symmetrical 
with respect to substitution q\ <-> — q 2 . It is quite a non-trivial task to notice a reason 
for this symmetry in the starting expression ( |103|) . Nevertheless, this symmetry is hidden 
there. It is a consequence of the invariance of the expression inside the curly brackets 
in ( |103|) , as well as the phase space element dx/(x(l — x))dli , under the "left-right 
symmetry" transformation ]T8| , i.e. under the substitution: 



-q 2 , h <-+l 2 , x <r+ 
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(107) 



The quark-antiquark pair production contribution was calculated in Refs. [20]. The 
result is 
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4 BFKL Pomeron in the next-to-leading approxima- 
tion 



Now we have all contributions to the NLLA kernel. Remind, that it has the same form 
fl52|) as the LLA one, where the gluon Regge trajectory in the two-loop approximation 
is given by ([H]) and the part K, r (qi,q 2 ) related with the real particle production (|65|) 
contains the contributions from the one-gluon (B3D, two-gluon (|105|) and quark-antiquark 



( |108| ) production in the Reggeon-Reggeon collisions. Representing this part as the sum 
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of the LLA contribution, expressed in terms of the renormalized coupling constant, and 
the one-loop correction JC^(q\, <f 2 ), we have for this correction fl2jj : 
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The remarkable fact exhibited by Eq. ( |110| ) is the cancellation of the infrared singularities 
[ Em (remind, that separate terms in (|B"5|) are singular as 1/e 2 ) at fixed fci = g*i — <f 2 , 
where we can expand (A; 2 //i 2 ) e in powers of e. This expansion is not performed in (|110[) 
because for the terms having singularity at k 2 = the region of such small fci, that is 
ln(l//c 1 2 ) ~ 1/e , does contribute in the integral over k±. The singular contributions given 
by this region are canceled, in turn, in the BFKL equation by the singular terms in the 
gluon trajectory 



As well as in LLA, only the kernel averaged over the angle between the momenta q\ 
and <f 2 is relevant until we don't consider spin correlations. For the averaged one-loop 
correction we have: 
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Instead of the dimensional regularization we can use the cut off |<fi 2— 9*2 2 | > ^ 2 changing 
correspondingly the virtual part. Using Eqs. (|52]), ([TCp, (|109 ) it is possible to verify that 
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the averaged NLLA kernel at e — > is equivalent to the expression 
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in the two-dimensional transverse space, with A — > 0. Of course, the dependence from A 
disappears when the kernel acts on a function. Moreover, the representation (|113|) permits 
to find such form of the kernel for which this cancellation is evident, just in the same way 
as it was done at transition from (f)6D to (1591). We obtain: 
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The [a - dependence in the right hand side of this equality leads to the violation of the 
scale invariance and is related with running the QCD coupling constant. 



The form ( |1 14| ) is very convenient for finding the action of the kernel on the eigen- 
functions q^ 1 of the Born kernel: 
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were within our accuracy we expressed the result in terms of the running coupling constant 
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For the relative correction r{^) defined as x^\l) = ~ r (l)x B (l) m the symmetrical 
point 7 = 1/2, corresponding to the eigenfunction of the LLA kernel with the largest 
eigenvalue, we have |21| 
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It shows that the correction is very large. 



6,46 + 0.05^ + 2.66^. (119) 



In some sense, the large value of the correction is natural and is a consequence of 
the large value of the LLA Pomeron intercept ujp = 4iVm2 a s (q 2 )/7[ . If we express the 
corrected intercept ujp in terms of the Born one, we obtain 
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(120) 



The coefficient 2.4 does not look very large. Moreover, it corresponds to the rapidity 
interval where correlations become important in the hadron production processes. 

Nevertheless, these numerical estimates show, that in the kinematical region of HERA 
probably it is not enough to take into account only the next-to-leading correction. For 
example, if a s (q 2 ) = 0.15, where the Born intercept is Up = 4Na s (q 2 )/ir ln2 = .39714, 
the relative correction for nf = is very big: 
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The maximal value of ujp ~ 0.1 is obtained for a s (q 2 ) ~ 0.08. 

But it is necessary to realize that, firstly, the above estimates are quite straightfor- 
ward and do not take into account neither the influence of the running coupling on the 
eigenfunctions nor the nonperturbative effects J23|; secondly, the value of the correction 
strongly depends on its representation. For example, if one takes into account the next- 
to-leading correction by the corresponding increase of the argument of the running QCD 
coupling constant, Up at a s (q 2 ) = 0.15 turns out to be only two times smaller, than its 
Born value. 

The above results can be applied for the calculation of anomalous dimensions of the 
local operators in the vicinity of the point uj = J — 1 = 0. The deep-inelastic moments 
JF9 B (g 2 ) defined as follows 

^V)=/JTv(^) "^(^s), (122) 

where 
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obey the integral equation of the type fl5~ll) with the inhomogeneous term equal to $s(g)/ (2nq 
and the kernel 
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The action of the modified kernel on the Born eigenfunctions q 2 2 ^' can t> e calculated 
easily: 
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where 

X {1) (7) = X (1) (7)-^(7)X'(7)- (126) 

The anomalous dimensions 7a, (a s ) = 7o(a s /u;) +a s 7i(a s /w) of the twist-2 operators near 
point uj = are determined from the solution of the equation 



„ = zdL ( X B (7) + ^ x -« (7) ) _ sS* (li±^! + 

7T \ 7T / 7T 7 471"^ \ 3 7^ 



%(10 + 13/iV 2 ) + 3£5 _ _ 117^ n, 71 _ _ 

9 7 iV + 27 WJ 3 6 + N* [ 27 9 } ) {U<) 

for 7 — > 0. For the low orders of the perturbation theory we reproduce the known results 
and predict the higher loop correction for uj —>■ 0: 

a s NA 11 tx f . / oi s \ 2 Tit N ,5 13, 



7T v cj 12 6iV 3 ' V 7T / 6cj v 3 6iV 2 



27 



2C(3) + . (128) 



4uo 2 V 7T / \ 27 ^ v ' 3 6 A 3 '27 9 



The results presented in this Section were obtained in Ref.|2T[. They differ from the 



corresponding results of Ref.p4f because the results of pi] were obtained for so called 



"irreducible part" of the kernel and therefore are incomplete. After appearance of Ref.pT 
the results obtained there were confirmed in Ref. P5fl . 

Finally, let us estimate the change of the maximal value of the anomalous dimension 
lw{oc 3 ) which is determined by the position of the extremum of the function 0^(7) ( 127 ). In 



LLA this extremum is reached at j^ax = 1/2 independently from a s . In NLLA, assuming 
the validity of the perturbative expansion, we obtain 

lmax ~ 2 tt ( X B a))" 2 + 4tt 1 6 3N + ^ j - 1 yj 
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